Vittrup Man-The life-history of a genetic foreigner in Neolithic Denmark.

The lethally maltreated body of Vittrup Man was deposited in a Danish bog, probably as part of a ritualised sacrifice. It happened between c. 3300 and 3100 cal years BC, i.e., during the period of the local farming-based Funnel Beaker Culture. In terms of skull morphological features, he differs from the majority of the contemporaneous farmers found in Denmark, and associates with hunter-gatherers, who inhabited Scandinavia during the previous millennia. His skeletal remains were selected for transdisciplinary analysis to reveal his life-history in terms of a population historical perspective. We report the combined results of an integrated set of genetic, isotopic, physical anthropological and archaeological analytical approaches. Strontium signature suggests a foreign birthplace that could be in Norway or Sweden. In addition, enamel oxygen isotope values indicate that as a child he lived in a colder climate, i.e., to the north of the regions inhabited by farmers. Genomic data in fact demonstrates that he is closely related to Mesolithic humans known from Norway and Sweden. Moreover, dietary stable isotope analyses on enamel and bone collagen demonstrate a fisher-hunter way of life in his childhood and a diet typical of farmers later on. Such a variable life-history is also reflected by proteomic analysis of hardened organic deposits on his teeth, indicating the consumption of forager food (seal, whale and marine fish) as well as farmer food (sheep/goat). From a dietary isotopic transect of one of his teeth it is shown that his transfer between societies of foragers and farmers took place near to the end of his teenage years.

Elevated genetic risk for multiple sclerosis emerged in steppe pastoralist populations.

Multiple sclerosis (MS) is a neuro-inflammatory and neurodegenerative disease that is most prevalent in Northern Europe. Although it is known that inherited risk for MS is located within or in close proximity to immune-related genes, it is unknown when, where and how this genetic risk originated1. Here, by using a large ancient genome dataset from the Mesolithic period to the Bronze Age2, along with new Medieval and post-Medieval genomes, we show that the genetic risk for MS rose among pastoralists from the Pontic steppe and was brought into Europe by the Yamnaya-related migration approximately 5,000 years ago. We further show that these MS-associated immunogenetic variants underwent positive selection both within the steppe population and later in Europe, probably driven by pathogenic challenges coinciding with changes in diet, lifestyle and population density. This study highlights the critical importance of the Neolithic period and Bronze Age as determinants of modern immune responses and their subsequent effect on the risk of developing MS in a changing environment.

100 ancient genomes show repeated population turnovers in Neolithic Denmark.

Major migration events in Holocene Eurasia have been characterized genetically at broad regional scales1-4. However, insights into the population dynamics in the contact zones are hampered by a lack of ancient genomic data sampled at high spatiotemporal resolution5-7. Here, to address this, we analysed shotgun-sequenced genomes from 100 skeletons spanning 7,300 years of the Mesolithic period, Neolithic period and Early Bronze Age in Denmark and integrated these with proxies for diet (13C and 15N content), mobility (87Sr/86Sr ratio) and vegetation cover (pollen). We observe that Danish Mesolithic individuals of the Maglemose, Kongemose and Ertebølle cultures form a distinct genetic cluster related to other Western European hunter-gatherers. Despite shifts in material culture they displayed genetic homogeneity from around 10,500 to 5,900 calibrated years before present, when Neolithic farmers with Anatolian-derived ancestry arrived. Although the Neolithic transition was delayed by more than a millennium relative to Central Europe, it was very abrupt and resulted in a population turnover with limited genetic contribution from local hunter-gatherers. The succeeding Neolithic population, associated with the Funnel Beaker culture, persisted for only about 1,000 years before immigrants with eastern Steppe-derived ancestry arrived. This second and equally rapid population replacement gave rise to the Single Grave culture with an ancestry profile more similar to present-day Danes. In our multiproxy dataset, these major demographic events are manifested as parallel shifts in genotype, phenotype, diet and land use.

The selection landscape and genetic legacy of ancient Eurasians.

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

Population genomics of post-glacial western Eurasia.

Western Eurasia witnessed several large-scale human migrations during the Holocene1-5. Here, to investigate the cross-continental effects of these migrations, we shotgun-sequenced 317 genomes-mainly from the Mesolithic and Neolithic periods-from across northern and western Eurasia. These were imputed alongside published data to obtain diploid genotypes from more than 1,600 ancient humans. Our analyses revealed a 'great divide' genomic boundary extending from the Black Sea to the Baltic. Mesolithic hunter-gatherers were highly genetically differentiated east and west of this zone, and the effect of the neolithization was equally disparate. Large-scale ancestry shifts occurred in the west as farming was introduced, including near-total replacement of hunter-gatherers in many areas, whereas no substantial ancestry shifts happened east of the zone during the same period. Similarly, relatedness decreased in the west from the Neolithic transition onwards, whereas, east of the Urals, relatedness remained high until around 4,000 BP, consistent with the persistence of localized groups of hunter-gatherers. The boundary dissolved when Yamnaya-related ancestry spread across western Eurasia around 5,000 BP, resulting in a second major turnover that reached most parts of Europe within a 1,000-year span. The genetic origin and fate of the Yamnaya have remained elusive, but we show that hunter-gatherers from the Middle Don region contributed ancestry to them. Yamnaya groups later admixed with individuals associated with the Globular Amphora culture before expanding into Europe. Similar turnovers occurred in western Siberia, where we report new genomic data from a 'Neolithic steppe' cline spanning the Siberian forest steppe to Lake Baikal. These prehistoric migrations had profound and lasting effects on the genetic diversity of Eurasian populations.

The genomic history of the indigenous people of the Canary Islands.

The indigenous population of the Canary Islands, which colonized the archipelago around the 3rd century CE, provides both a window into the past of North Africa and a unique model to explore the effects of insularity. We generate genome-wide data from 40 individuals from the seven islands, dated between the 3rd-16rd centuries CE. Along with components already present in Moroccan Neolithic populations, the Canarian natives show signatures related to Bronze Age expansions in Eurasia and trans-Saharan migrations. The lack of gene flow between islands and constant or decreasing effective population sizes suggest that populations were isolated. While some island populations maintained relatively high genetic diversity, with the only detected bottleneck coinciding with the colonization time, other islands with fewer natural resources show the effects of insularity and isolation. Finally, consistent genetic differentiation between eastern and western islands points to a more complex colonization process than previously thought.

Imputation of ancient human genomes.

Due to postmortem DNA degradation and microbial colonization, most ancient genomes have low depth of coverage, hindering genotype calling. Genotype imputation can improve genotyping accuracy for low-coverage genomes. However, it is unknown how accurate ancient DNA imputation is and whether imputation introduces bias to downstream analyses. Here we re-sequence an ancient trio (mother, father, son) and downsample and impute a total of 43 ancient genomes, including 42 high-coverage (above 10x) genomes. We assess imputation accuracy across ancestries, time, depth of coverage, and sequencing technology. We find that ancient and modern DNA imputation accuracies are comparable. When downsampled at 1x, 36 of the 42 genomes are imputed with low error rates (below 5%) while African genomes have higher error rates. We validate imputation and phasing results using the ancient trio data and an orthogonal approach based on Mendel's rules of inheritance. We further compare the downstream analysis results between imputed and high-coverage genomes, notably principal component analysis, genetic clustering, and runs of homozygosity, observing similar results starting from 0.5x coverage, except for the African genomes. These results suggest that, for most populations and depths of coverage as low as 0.5x, imputation is a reliable method that can improve ancient DNA studies.

The transcriptional and regulatory identity of erythropoietin producing cells.

Erythropoietin (Epo) is the master regulator of erythropoiesis and oxygen homeostasis. Despite its physiological importance, the molecular and genomic contexts of the cells responsible for renal Epo production remain unclear, limiting more-effective therapies for anemia. Here, we performed single-cell RNA and transposase-accessible chromatin (ATAC) sequencing of an Epo reporter mouse to molecularly identify Epo-producing cells under hypoxic conditions. Our data indicate that a distinct population of kidney stroma, which we term Norn cells, is the major source of endocrine Epo production in mice. We use these datasets to identify the markers, signaling pathways and transcriptional circuits characteristic of Norn cells. Using single-cell RNA sequencing and RNA in situ hybridization in human kidney tissues, we further provide evidence that this cell population is conserved in humans. These preliminary findings open new avenues to functionally dissect EPO gene regulation in health and disease and may serve as groundwork to improve erythropoiesis-stimulating therapies.

A 2-million-year-old ecosystem in Greenland uncovered by environmental DNA.

Late Pliocene and Early Pleistocene epochs 3.6 to 0.8 million years ago1 had climates resembling those forecasted under future warming2. Palaeoclimatic records show strong polar amplification with mean annual temperatures of 11-19 °C above contemporary values3,4. The biological communities inhabiting the Arctic during this time remain poorly known because fossils are rare5. Here we report an ancient environmental DNA6 (eDNA) record describing the rich plant and animal assemblages of the Kap København Formation in North Greenland, dated to around two million years ago. The record shows an open boreal forest ecosystem with mixed vegetation of poplar, birch and thuja trees, as well as a variety of Arctic and boreal shrubs and herbs, many of which had not previously been detected at the site from macrofossil and pollen records. The DNA record confirms the presence of hare and mitochondrial DNA from animals including mastodons, reindeer, rodents and geese, all ancestral to their present-day and late Pleistocene relatives. The presence of marine species including horseshoe crab and green algae support a warmer climate than today. The reconstructed ecosystem has no modern analogue. The survival of such ancient eDNA probably relates to its binding to mineral surfaces. Our findings open new areas of genetic research, demonstrating that it is possible to track the ecology and evolution of biological communities from two million years ago using ancient eDNA.

Genomic ancestry, diet and microbiomes of Upper Palaeolithic hunter-gatherers from San Teodoro cave.

Recent improvements in the analysis of ancient biomolecules from human remains and associated dental calculus have provided new insights into the prehistoric diet and genetic diversity of our species. Here we present a multi-omics study, integrating metagenomic and proteomic analyses of dental calculus, and human ancient DNA analysis of the petrous bones of two post-Last Glacial Maximum (LGM) individuals from San Teodoro cave (Italy), to reconstruct their lifestyle and the post-LGM resettlement of Europe. Our analyses show genetic homogeneity in Sicily during the Palaeolithic, representing a hitherto unknown Italian genetic lineage within the previously identified Villabruna cluster. We argue that this lineage took refuge in Italy during the LGM, followed by a subsequent spread to central-western Europe. Analysis of dental calculus showed a diet rich in animal proteins which is also reflected on the oral microbiome composition. Our results demonstrate the power of this approach in the study of prehistoric humans and will enable future research to reach a more holistic understanding of the population dynamics and ecology.

A Middle Pleistocene Denisovan molar from the Annamite Chain of northern Laos.

The Pleistocene presence of the genus Homo in continental Southeast Asia is primarily evidenced by a sparse stone tool record and rare human remains. Here we report a Middle Pleistocene hominin specimen from Laos, with the discovery of a molar from the Tam Ngu Hao 2 (Cobra Cave) limestone cave in the Annamite Mountains. The age of the fossil-bearing breccia ranges between 164-131 kyr, based on the Bayesian modelling of luminescence dating of the sedimentary matrix from which it was recovered, U-series dating of an overlying flowstone, and U-series-ESR dating of associated faunal teeth. Analyses of the internal structure of the molar in tandem with palaeoproteomic analyses of the enamel indicate that the tooth derives from a young, likely female, Homo individual. The close morphological affinities with the Xiahe specimen from China indicate that they belong to the same taxon and that Tam Ngu Hao 2 most likely represents a Denisovan.

Ancient Human Genomes and Environmental DNA from the Cement Attaching 2,000-Year-Old Head Lice Nits.

Over the past few decades, there has been a growing demand for genome analysis of ancient human remains. Destructive sampling is increasingly difficult to obtain for ethical reasons, and standard methods of breaking the skull to access the petrous bone or sampling remaining teeth are often forbidden for curatorial reasons. However, most ancient humans carried head lice and their eggs abound in historical hair specimens. Here we show that host DNA is protected by the cement that glues head lice nits to the hair of ancient Argentinian mummies, 1,500-2,000 years old. The genetic affinities deciphered from genome-wide analyses of this DNA inform that this population migrated from north-west Amazonia to the Andes of central-west Argentina; a result confirmed using the mitochondria of the host lice. The cement preserves ancient environmental DNA of the skin, including the earliest recorded case of Merkel cell polyomavirus. We found that the percentage of human DNA obtained from nit cement equals human DNA obtained from the tooth, yield 2-fold compared with a petrous bone, and 4-fold to a bloodmeal of adult lice a millennium younger. In metric studies of sheaths, the length of the cement negatively correlates with the age of the specimens, whereas hair linear distance between nit and scalp informs about the environmental conditions at the time before death. Ectoparasitic lice sheaths can offer an alternative, nondestructive source of high-quality ancient DNA from a variety of host taxa where bones and teeth are not available and reveal complementary details of their history.

Transcriptional Activity of Predominant Streptococcus Species at Multiple Oral Sites Associate With Periodontal Status.

Background: Streptococcus species are predominant members of the oral microbiota in both health and diseased conditions. The purpose of the present study was to explore if different ecological characteristics, such as oxygen availability and presence of periodontitis, associates with transcriptional activity of predominant members of genus Streptococcus. We tested the hypothesis that genetically closely related Streptococcus species express different transcriptional activities in samples collected from environments with critically different ecological conditions determined by site and inflammatory status. Methods: Metagenomic and metatranscriptomic data was retrieved from 66 oral samples, subgingival plaque (n=22), tongue scrapings (n=22) and stimulated saliva (n=22) collected from patients with periodontitis (n=11) and orally healthy individuals (n=11). Species-specific transcriptional activity was computed as Log2(RNA/DNA), and transcriptional activity of predominant Streptococcus species was compared between multiple samples collected from different sites in the same individual, and between individuals with different oral health status. Results: The predominant Streptococcus species were identified with a site-specific colonization pattern of the tongue and the subgingival plaque. A total of 11, 4 and 2 pathways expressed by S. parasanguinis, S. infantis and S. salivarius, respectively, were recorded with significantly higher transcriptional activity in saliva than in tongue biofilm in healthy individuals. In addition, 18 pathways, including pathways involved in synthesis of peptidoglycan, amino acid biosynthesis, glycolysis and purine nucleotide biosynthesis expressed by S. parasanguinis and 3 pathways expressed by S. salivarius were identified with significantly less transcriptional activity in patients with periodontitis. Conclusion: Data from the present study significantly demonstrates the association of site-specific ecological conditions and presence of periodontitis with transcriptional activity of the predominant Streptococcus species of the oral microbiota. In particular, pathways expressed by S. parasanguinis being involved in peptidoglycan, amino acid biosynthesis, glycolysis, and purine nucleotide biosynthesis were identified to be significantly associated with oral site and/or inflammation status.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears.

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

Genomic Steppe ancestry in skeletons from the Neolithic Single Grave Culture in Denmark.

The Gjerrild burial provides the largest and best-preserved assemblage of human skeletal material presently known from the Single Grave Culture (SGC) in Denmark. For generations it has been debated among archaeologists if the appearance of this archaeological complex represents a continuation of the previous Neolithic communities, or was facilitated by incoming migrants. We sampled and analysed five skeletons from the Gjerrild cist, buried over a period of c. 300 years, 2600/2500-2200 cal BCE. Despite poor DNA preservation, we managed to sequence the genome (>1X) of one individual and the partial genomes (0.007X and 0.02X) of another two individuals. Our genetic data document a female (Gjerrild 1) and two males (Gjerrild 5 + 8), harbouring typical Neolithic K2a and HV0 mtDNA haplogroups, but also a rare basal variant of the R1b1 Y-chromosomal haplogroup. Genome-wide analyses demonstrate that these people had a significant Yamnaya-derived (i.e. steppe) ancestry component and a close genetic resemblance to the Corded Ware (and related) groups that were present in large parts of Northern and Central Europe at the time. Assuming that the Gjerrild skeletons are genetically representative of the population of the SGC in broader terms, the transition from the local Neolithic Funnel Beaker Culture (TRB) to SGC is not characterized by demographic continuity. Rather, the emergence of SGC in Denmark was part of the Late Neolithic and Early Bronze Age population expansion that swept across the European continent in the 3rd millennium BCE, resulting in various degrees of genetic replacement and admixture processes with previous Neolithic populations.

Population genomics of the Viking world.

The maritime expansion of Scandinavian populations during the Viking Age (about AD 750-1050) was a far-flung transformation in world history1,2. Here we sequenced the genomes of 442 humans from archaeological sites across Europe and Greenland (to a median depth of about 1×) to understand the global influence of this expansion. We find the Viking period involved gene flow into Scandinavia from the south and east. We observe genetic structure within Scandinavia, with diversity hotspots in the south and restricted gene flow within Scandinavia. We find evidence for a major influx of Danish ancestry into England; a Swedish influx into the Baltic; and Norwegian influx into Ireland, Iceland and Greenland. Additionally, we see substantial ancestry from elsewhere in Europe entering Scandinavia during the Viking Age. Our ancient DNA analysis also revealed that a Viking expedition included close family members. By comparing with modern populations, we find that pigmentation-associated loci have undergone strong population differentiation during the past millennium, and trace positively selected loci-including the lactase-persistence allele of LCT and alleles of ANKA that are associated with the immune response-in detail. We conclude that the Viking diaspora was characterized by substantial transregional engagement: distinct populations influenced the genomic makeup of different regions of Europe, and Scandinavia experienced increased contact with the rest of the continent.

Ancient Jomon genome sequence analysis sheds light on migration patterns of early East Asian populations.

Anatomically modern humans reached East Asia more than 40,000 years ago. However, key questions still remain unanswered with regard to the route(s) and the number of wave(s) in the dispersal into East Eurasia. Ancient genomes at the edge of the region may elucidate a more detailed picture of the peopling of East Eurasia. Here, we analyze the whole-genome sequence of a 2,500-year-old individual (IK002) from the main-island of Japan that is  characterized with a typical Jomon culture. The phylogenetic analyses support multiple waves of migration, with IK002 forming a basal lineage to the East and Northeast Asian genomes examined, likely representing some of the earliest-wave migrants who went north from Southeast Asia to East Asia. Furthermore, IK002 shows strong genetic affinity with the indigenous Taiwan aborigines, which may support a coastal route of the Jomon-ancestry migration. This study highlights the power of ancient genomics to provide new insights into the complex history of human migration into East Eurasia.

Diverse variola virus (smallpox) strains were widespread in northern Europe in the Viking Age.

Smallpox, one of the most devastating human diseases, killed between 300 million and 500 million people in the 20th century alone. We recovered viral sequences from 13 northern European individuals, including 11 dated to ~600-1050 CE, overlapping the Viking Age, and reconstructed near-complete variola virus genomes for four of them. The samples predate the earliest confirmed smallpox cases by ~1000 years, and the sequences reveal a now-extinct sister clade of the modern variola viruses that were in circulation before the eradication of smallpox. We date the most recent common ancestor of variola virus to ~1700 years ago. Distinct patterns of gene inactivation in the four near-complete sequences show that different evolutionary paths of genotypic host adaptation resulted in variola viruses that circulated widely among humans.